Constraint on the Post-Newtonian Parameter 7 on Galactic Size Scales 
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We constrain the post-Newtonian gravity parameter 7 on kiloparsec scales by comparing the 
masses of 15 eUiptical lensing galaxies from the Sloan Lens ACS Survey as determined in two inde- 
pendent ways. The first method assumes only that Newtonian gravity is correct and is independent 
of 7, while the second uses gravitational lensing which depends on 7. More specifically, we combine 
Einstein radii and radial surface-brightness gradient measurements of the lens galaxies with em- 
pirical distributions for the mass concentration and velocity anisotropy of elliptical galaxies in the 
local universe to predict 7-dependent probability distributions for the lens-galaxy velocity disper- 
sions. By comparing with observed velocity dispersions, we derive a maximum-likelihood value of 
7 = 0.98 ±0.07 (68% confidence). This result is in excellent agreement with the prediction of general 
relativity that has previously been verified to this accuracy only on solar-system length scales. 
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Gravitational lenses provide some of the most inspi- 
rational and thought-provoking images in astronomy. R 
is often said in the popular press that these represent a 
wonderful verification of Einstein's theory of general rela- 
tivity (GR). In 1937, Zwicky 1] proposed that lensing by 
distant galaxies and clusters would furnish both a test of 
GR and a tool for the measurement of the lensing masses. 
Since the discovery of the first gravitational lenses, how- 
ever, astronomers have emphasized the latter application 
over the former. In fact, gravitational lensing alone can- 
not separately determine the masses of lenses and test 
the weak-field limit of the Schwarzschild metric which 
underlies the theory of lensing. 

Since Eddington's original solar eclipse expedition of 
1919, the Schwarzschild metric has been extensively 
probed in the weak field limit within the solar sys- 
tem (e.g. j2]) and with binary radio pulsars (e.g. Q])- 
In all cases, however, the scales involved are of order 
light seconds. Most recently, the Cassini mission has de- 
termined the post-Newtonian parameter 7 (e.g. 0, 0), 
described below, by directly measuring the Shapiro de- 
lay as radio signals pass by the Sun in their travel 
from the spacecraft to the Earth, giving a measured 
value of 7 = 1 (2.1 ± 2.3) x IQ-^ 0. In this paper 
we set precise constraints on 7 on galactic scale sizes 
of several kiloparsecs (1 parsec ~ 3.1 x 10^^ cm): im- 
pact parameters ~ 10^^ times larger than in Sun-grazing 
solar-system tests. We accomplish this by comparing the 
masses of relatively low-redshift (0.06 < z < 0.33) galax- 
ies as deduced from strong gravitational lensing with 
their masses as estimated from their stellar orbital dy- 
namics. Since the former determination depends upon 
GR through the value of 7, whereas the latter depends 
only upon Newtonian gravity, the quantitative agreement 
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between the two methods can place a direct constraint 
upon 7. Our analysis constrains 7 to have a value of 
0.98 ±0.07 (68% confidence), in excellent agreement with 
the GR prediction of 7 = 1. This statistical precision 
is enabled by data for a large and homogeneous sample 
of recently discovered gravitational lens galaxies. Simi- 
lar analyses have been carried out previously, but with 
much lower statistical significance. Nottale [3 applied 
this type of test to the lensing galaxy cluster A370 to ob- 
tain 0.87 < 7 < 1.55, subject to assumptions about the 
mass structure of the cluster. Dar Q found agreement 
within ^30% error bars between the observed and pre- 
dicted image separations in a heterogeneous sample of 5 
gravitational-lens systems, also for fixed mass-structure 
assumptions, but did not translate this result into a con- 
straint on 7. Sirousse-Zia 01 developed formalism for 
the test, expressing the dependence upon 7 of the lens- 
ing properties of the singular isothermal sphere galaxy 
model (see below), but did not derive any observational 
results. 

The dependence of gravitational lensing upon 7 can 
be derived from a general form of the metric for a point 
mass m in the weak field limit for a space-time that yields 
Newtonian gravity: 

dr^ = dt' (1 - V) - (1 + ^) - ^"dc^" ■ (1) 

The parameter 7 equals unity for the Schwarzschild met- 
ric, (f) is the angle in the invariant orbital plane, and G 
and c have been set equal to 1. In this formulation, the 
weak-field gravitational acceleration can be calculated, 
and is verified to be — Gm/r^ as given by Newton and 
independent of the parameter 7. From the speed of light 
in the radial direction as inferred by an external ob- 
server at infinity, one may define an effective index of 
refraction for the space surrounding the point mass, and 
hence around any arbitrary mass distribution through 
the principle of superposition. This quantity is used to 
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compute the Shapiro delay, and for sufficiently isolated 
gravitational-lens systems one may invoke the thin-lens 
approximation. Thus gravitational lensing can be for- 
mulated in terms of two-dimensional Fermat time delay 
surfaces ,10. 11] . By extremizing the sum of the Shapiro 
and geometric time delays, one obtains the "lens equa- 
tion" that relates position in the observed image plane 
to location in the unobserved and unlensed source plane: 

es = e-^^-^wm ■ (2) 

Here, 6s is the angular source location, 6 is the angular 
location of the image, and is a scaled line-of-sight 
integral of the Newtonian gravitational potential of the 
lensing object (see Eq. 48 of for the explicit defini- 
tion of i} in this context). The Einstein radius, defined 
by 0s = for circularly-symmetric projected potentials, 
simply scales with the factor (1 -|-7)/2. The difficulty in 
inferring anything about the 7 parameter is that it ap- 
pears only as a product with the relevant masses of the 
problem. Unless one knows the lensing mass via some 
method other than lensing, there is no distinction be- 
tween, e.g., the solution for a lensing mass M and 7=1 
and the solution for a lensing mass of 2M and 7 = 0. 

The actual observational test that we perform here 
can be understood most simply in terms of the singu- 
lar isothermal sphere (SIS) galaxy model, though we in 
fact allow for more general galaxy models in our anal- 
ysis. The SIS is a dynamically self-consistent, spher- 
ically symmetric three-dimensional profile with density 
p{r) = a1 / {2-KGr'^) . It is perhaps the best one-parameter 
model for elliptical galaxies, characterized by an isotropic 
and radially constant dispersion of stellar orbital veloci- 
ties CTu that determines the mass within any given radius. 
Through gravitational lensing, the SIS forms a ring image 
of any background objects along the same line of sight, 
with an angular "Einstein radius" of 

eE = {l+i)2^{<yl/^){DLs/Ds) , (3) 

where D^s a-nd Ds are distance measures discussed fur- 
ther below. Thus, within the context of the SIS model, 
measurements of the Einstein radii and velocity disper- 
sions of gravitational lens galaxies can be used to con- 
strain the 7 parameter . (From an observational stand- 
point, velocity dispersion is defined in this work as the 
luminosity- weighted integral of the second moment of the 
stellar velocity distribution along the observational line of 
sight. It may be measured from the spectrum of a galaxy 
whose stellar population is spatially unresolved by fitting 
for the broadening of stellar atomic absorption lines that 
best reproduces the features seen in the integrated galaxy 
spectrum.) 

Spatially resolved, high signal-to-noise observations of 
the kinematics of elliptical galaxies, whose stellar popu- 
lations are dynamically "hot" (i.e., not characterized by 
ordered circular orbits) , can be used in combination with 
detailed modeling to deduce radial density profiles and to 



directly test the validit y of the isothermal approximation. 
Such studies (e.g. El El El El El El 1 9] ) indicate a 
modest amount of dark matter within one half-light ra- 
dius and an approximately isothermal density profile in 
elliptical galaxies that are nearby enough to permit the 
necessary observations. Gravitational lens galaxies (here- 
after GLGs), by contrast, are typically elliptical galaxies 
at relatively high redshift and with significant contami- 
nation from the light of lensed quasar images, and thus 
do not allow the same type of dynamical analysis as local 
ellipticals. Thus lensing itself is usually the only robust 
measurement of the mass of the GLG, and a model for 
the density profile must simply be assumed. An excep- 
tion is the Lenses Structure and Dynamics [LSD] Survey 
El m m m m, which has used lens-galaxy velocity 
dispersions in combination with Einstein radii and the 
implicit GR assumption of 7 = 1 to place constraints on 
the relative contributions of luminous and dark matter 
in distant GLGs. 

The current work is enabled by the significant new 
sample of relatively low redshift (0.06 < z < 0.33) ellip- 
tical GLGs presented by the Sloan Lens ACS (SLAGS) 
Survey jH M El- The SLAGS GLGs, which were 
discovered within the Sloan Digital Sky Survey spec- 
troscopic database (SDSS: 28]), are particularly distin- 
guished by the ease and accuracy with which their sur- 
face brightness profiles and stellar dynamics can be mea- 
sured, relative to previously known GLGs. Publicly avail- 
able SDSS velocity-dispersion measurements of SLAGS 
GLGs provide a single, lensing-independent dynamical 
mass scale for each system, though they do not allow a 
lensing-independent determination of the GLG density 
profiles. In this work we adopt the hypothesis that, by 
virtue of their relatively low redshifts, SLAGS GLGs are 
sufficiently like nearby elliptical galaxies for their dis- 
tribution in density profiles to be approximated by the 
distribution deduced for nearby elliptical galaxies from 
the application of detailed dynamical modeling. The 
masses of SLAGS GLGs can thus be inferred indepen- 
dently of lensing from their SDSS velocity dispersions, 
with an uncertainty quantified by both measurement er- 
ror and intrinsic scatter in density profile (and velocity 
anisotropy, discussed below) . The most significant evolu- 
tion to be expected between nearby elliptical galaxies and 
the SLAGS GLGs is in luminosity, which we circumvent 
by working entirely with mass, shape, and dynamical ob- 
servablcs. 

In order to work with directly observable quantities (of 
which galaxy mass is not one), we frame our analysis in 
terms of a comparison between the stellar velocity disper- 
sions of GLGs (i) as observed and (ii) as predicted from 
their lensing Einstein radii for a given value of 7. Before 
proceeding into more detail, we note that the observed 
stellar dispersions of SLAGS GLGs agree within observa- 
tional errors with their Einstein radii when the latter are 
directly translated into velocity dispersions using Eq. O 
and assuming 7=1 (see Fig. 5 of El]). To incorpo- 
rate the possible effects of non-isothermal lens profiles. 
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we employ the self-similar, axisymmetric model described 
by Koopmans j29l |. This model approximates the mass- 
and luminosity-density profiles of GLGs with power-law 
forms: i.e. three-dimensional mass density p(r) oc r~" 
and luminosity density ^{r) cx (here we have used a 
to replace the 7' of |29] to avoid confusion with the post- 
Newtonian parameter of interest). With a and 5 con- 
sidered as separate parameters, the model corresponds 
to a scale-free galaxy with a constant logarithmic ra- 
dial mass-to-light ratio gradient, capable of describing 
an increasing dark-matter fraction at increasing radius. 
Larger values of a and S correspond to more centrally 
concentrated mass and light profiles. The observed an- 
gular Einstein radius of a GLG gives a measurement 
of the enclosed mass. Taking the mass normalization 
set by 9e and fixed values for a and S, the spherically- 
symmetric steady-state Jeans equation is integrated an- 
alytically for velocity dispersion in the radial direction 
as a function of radius. (The Jeans equation is obtained 
from the velocity moments of the coUisionless Boltzmann 
equation for the phase space distribution of stellar orbits 
in Newtonian gravity; see, e.g., [s^l) The model further 
assumes a constant radial profile for the conventionally 
defined velocity anisotropy parameter /3 which relates ve- 
locity dispersions in the radial and tangential directions: 
/9 = 1 — Cu tan/'^u rad- Jeans-cquatiou solution as- 

sumes the validity of Newtonian gravitation on the rele- 
vant scales, but is independent of 7 in the weak-field limit 
as discussed above. The luminosity-weighted squared ve- 
locity dispersion is integrated along the line of sight to 
give an analytic prediction for the observed velocity dis- 
persion as a function of position within the image of the 
GLG. The simplest SIS model is described by the special 
case of a = ^ = 2 and /? = 0. 

The sample we analyze consists of the 15 SLAGS GLGs 
with published Einstein radii, determined by fitting sin- 
gular isothermal ellipsoid (SIE) lens models with the 
normalization of 31] to Hubble Space Telescope {HST) 
imaging data ^Tlj. Though the SIE model assumes a 
particular radial mass profile, measured Einstein radii 
depend only weakly upon the details of the radial mass 
distribution of the GLG. We obtain i5 values by comput- 
ing average logarithmic surface-brightness profile slopes 
for the 15 systems through non-linear least-squares fitting 
of elliptical power-law luminosity models convolved with 
the instrumental point-spread function to the HST imag- 
ing data within a radius of 1.8 arcseconds about the cen- 
ter of each GLG. (A three-dimensional brightness profile 
j/(r) oc will have a projected two-dimensional profile 
I{R) oc R-^+\) Both the mass and light models fit- 
ted to the HST data include a projected major-to-minor 
axis-ratio parameter q to allow for ellipticity. To con- 
nect these models to the axisymmetric approximation of 
the analytic Jeans equation-based framework, we use the 
interchange 



iso-density or iso-brightness contour. Stellar velocity dis- 
persion values for each GLG are taken from the output 
of the Princeton Id spectroscopic pipeline analysis of the 
SDSS database ^3^. The SDSS spectrograph integrates 
all atmospherically blurred galaxy light within a circular 
optical- fiber aperture of radius 1.5 arcseconds on the sky. 
The fractional systematic velocity-dispersion inaccuracy 
incurred in a circularly-symmetric Jeans analysis of an 
oblate galaxy is calculated by [s^ to be small provided 
that one uses angularly-averaged velocity dispersions and 
Einstein radii, as we do in our calculations here. We add 
this error in quadrature to the observational velocity- 
dispersion error estimates for each lens, though it is small 
(median contribution of 2%) and has a negligible effect 
on our result. 

To determine suitable probability distributions for the 
mass concentration a and velocity anisotropy f3 of the 
SLAGS GLGs, we make use of the results of Gerhard et 
al. .15] ■ who compute the density structure for a sam- 
ple of nearby elliptical galaxies using dynamical observa- 
tions. From Fig. 1 of that work, we determine the change 
in circular velocity (or, equivalently, mass enclosed) be- 
tween 0.2 and 0.6 half-light radii for each of the 17 galax- 
ies with data over that range. The chosen radial range 
corresponds roughly to the range probed by the SLAGS 
GLG sample. We assign to each of these galaxies a log- 
arithmic radial density slope a that will give the same 
relative change in circular velocity. This gives a sample 
mean log-slope of (a) — 1.93 and an intrinsic RMS vari- 
ation of CTq — 0.08. From Fig. 5 of the same work, we 
calculate the average and RMS variation of the velocity 
anisotropy parameter /? (e.g. |3) to be (/?) = 0.18 and 
= 0.13 (excluding the outlier galaxies NGC 4636 and 
NGC 4486B). We model the probability distributions of 
a and /3 for the SLAGS GLGs as uncorrelated Gaussians 
with these parameters. Neither of these properties ex- 
hibit significant correlation with galaxy mass in . 

For comparison with the observed velocity disper- 
sions, we integrate the position-dependent velocity dis- 
persion predicted by lensing over the SDSS spectroscopic 
optical-fiber aperture. The combined effect of the circu- 
lar SDSS aperture (radius 1.5 arcseconds) and the typ- 
ical image quality of FWHM 2 arcseconds (due to at- 
mospheric blurring) is well approximated by a circular 
Gaussian weighted aperture with a FWHM of 2.8 arc- 
seconds, which we combine with luminosity weighting in 
an integration of Eq. 2.4 of We furthermore incor- 
porate the dependence of gravitational lensing upon 7 
and eliminate the explicit dependence on the mass of the 
GLG in favor of the Einstein radius to obtain 
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which conserves the total mass or light within a given 



Here, the first occurrence of 9e is to be expressed in 
radians. The dimensionless function f (a, 6,(3) is given 
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in [23. The quantity 9a is the "Gaussian sigma" of the 
spatial weighting aperture: i.e., 2.8 arcseconds divided by 
2.355. Factors on the second line arise from our use of a 
Gaussian integration aperture. Ds and -Dls are angular- 
diameter distances from the observer to the lensed galaxy 
and from the GLG to the lensed galaxy respectively, spec- 
ifying the angle subtended by a given physical distance 
transverse to the line of sight. These distances are func- 
tions of the redshifts of the two galaxies, which have 
been measured to high accuracy jfor all SLAGS systems 
using SDSS spectroscopic data 25j. Dg and Dj^s also 
depend upon the adopted cosmology, and in principle 
would be altered through the Robertson- Walker metric 
by any change in 7 (as noted by Q). However, such 
changes would necessarily be accompanied by changes 
in the matter-energy density of the universe in order to 
reproduce the shape of the luminosity distance-redshift 
relation empirically constrained by type la supernova ob- 
servations out to redshift z k, 1 ,3^ 35, 36, 37] . Thus al- 
though we shall compute Ds and D^s under the assump- 
tion of the currently favored cosmology with density pa- 
rameters (r^M, ^^a) = (0.3, 0.7) [sll, we regard this recipe 
as a suitable proxy for the empirical distance-redshift re- 
lation. 

We use Eq. [S] to determine numerically the probabil- 
ity distribution of the "true" aperture velocity dispersion 
(^fo s )^^^; which we denote below by a for convenience, 
for each GLG as a function of 7. This is done by com- 
puting a over a grid of all physically allowed a and (3 
values (i.e., those that give finite central mass and pre- 
dict non-negative cr^) at fixed {S,9e,0a), assigning dif- 
ferential weights over the grid as given by the assumed 
probabilities of a and /3, and sorting in a to determine the 
cumulative probability distribution. This computation is 
repeated over a grid of 7 values from 0-4 to derive the 
desired probability density Pi{(t\j) for each lens. The 
probability density for obtaining a particular measured 
SDSS velocity dispersion usdss at fixed 7 is then given 
by 

-p3(a'sDSs|7) = J da P2{asDSsW)Pi{'^\l) ■ (6) 

For P2(o'sDSs|o'), we assume a Gaussian distribution 
about the true value, with a width given by the mea- 
surement error for each system. Statistical errors on the 
measured 6 and 9e values are neglected, since they have 
an exceedingly small effect compared to the errors on the 
fSDSS values and the assumed intrinsic variation in a and 
(3. 

For each GLG i, given its measured ctsdss, Eq. |Hlmay 
be interpreted as a likelihood distribution for the param- 
eter 7. Since we expect the true value of 7 to be universal 
and thus the same for all GLGs in the sample, we can 
express its likelihood distribution given all the data and 
the prior assumptions as 

-^(7)= n P3\4Lsh) ■ (7) 

lenses i 




Post-Newtonion y 



FIG. 1: Normalized likelihood densities for the post- 
Newtonian parameter 7 derived from individual SLACS 
GLGs (grey), with joint likelihood for the entire sample 
(black). The joint distribution is approximately Gaussian, 
giving 7 = 0.98 ± 0.07 (68% confidence). 

Fig. n shows the individual factors in this likelihood 
product, as well as the joint normalized likelihood den- 
sity. Though the distribution of individual maximum- 
likelihood values seen in Fig. ^ appears slightly bimodal 
between 7 values greater than 1 and less than 1, this effect 
is not statistically significant: a Kolmogorov-Smirnov 
test of the distribution of individual 7 values indicates a 
reasonable 33% probability of being drawn from a Gaus- 
sian distribution about 7=1 with a width equal to the 
mean RMS width of the individual 7 distributions (ap- 
proximately 0.3). Furthermore, the individual 7 values 
are not significantly correlated with any lens-galaxy ob- 
servables. The joint distribution £(7) is described well 
by a Gaussian in 7, from which we derive a maximum- 
likelihood constraint of 7 = 0.98 ±0.07 (68% confidence). 

Our maximum-likelihood value for 7 is in excellent 
agreement with the GR-predicted value of 7 = 1. This 
result is an important quantitative test of the theory of 
gravitation on scales much larger than have been probed 
previously with solar-system experiments. The current 
analysis is enabled by the new sample of GLGs from the 
SLAGS Survey for which lens masses can (implicitly) be 
estimated through a combination of observables and rea- 
sonable assumptions independent of lensing itself. We 
note that a spatially resolved and/or higher signal-to- 
noise observational determination of the stellar dynamics 
of the SLAGS GLGs would remove some of the reliance 
on prior assumptions based on the local universe, thus 
making the result more robust. 
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